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We present a scheme for generating a synthetic magnetic field and spin-orbit coupling via Raman 
coupling in highly magnetic lanthanide atoms such as dysprosium. Employing these atoms offer 
several advantages for realizing strongly correlated states and exotic spinor phases. The large spin 
and narrow optical transitions of these atoms allow the generation of synthetic magnetic fields an 
order of magnitude larger than those in the alkalis, but with considerable reduction of the heating 
rate for equal Raman coupling. The effective hamiltonian of these systems differs from that of the 
alkalis' by an additional nematic coupling term, which leads to a phase transition in the dressed 
states as detuning varies. For high-spin condensates, spin-orbit coupling leads to a spatially periodic 
structure, which is described in Majorana representation by a set of points moving periodically on 
a unit sphere. We name this a "Majorana spinor helix" in analogy to the persistent spin-| helix 
observed in electronic systems. 



In the past few years, several groups have realized a 
synthetic magnetic field either in traps or in optical lat- 
tices [IHl], and spin-orbit (SO) coupling P^TT] with al- 
kali atoms. These developments have highlighted intrigu- 
ing physics in the ultracold atomic gas context. Vor- 
tices and the classical Hall effect have been observed 
with a Bose condensate exposed to a synthetic magnetic 
field [11 [2]. SO-coupling in a Bose gas can lead to super- 
fluid phases with stripe order [12] [13] and a rich phase 
diagram [5[ [131 [H]? as well as modify the effective in- 
teraction between dressed-state atoms [6 . In addition, 
SO-coupling also leads to a divergent spin susceptibility, 
and the magnetic transition it implies has recently been 
observed [8l|T5]. Recently, SO-coupled Fermi gases are 
found to display interesting spin dynamics, topological 
transitions of Fermi surfaces, and spin-dependent band 
structure piTol. 

However, there are serious challenges with creating ex- 
otic quantum matter using the current scheme of generat- 
ing synthetic gauge fields. The small fine-structure split- 
ting of the excited level used in the Raman-coupling pro- 
cess for alkalis leads to significant heating through spon- 
taneous emission [16 . While lowering the laser intensity 
will reduce heating, it will also reduce the strength of 
the synthetic gauge field, pushing the high-field regime 
of novel correlated physics beyond reach. As of now, the 
number of vortices generated by synthetic magnetic field 
is far below that generated by rotations [l]. 

We suggest the use of lanthanide atoms such as 
Dy [TTl [18] and Er \19\ to overcome these challenges. 
The particular atomic structure of these atoms — narrow 
linewidth transitions, large ground-state orbital and spin 
angular momenta, and large fine-structure splitting — 
offer many advantages over the alkalis. As we explain 
later, the narrow-line transitions (2-kHz wide in Dy ver- 
sus 6 MHz in Rb) and the L > nature of the Dy ground 



and excited states provide a significant increase in Raman 
coupling without additional heating, or conversely, much 
less heating at a fixed Raman coupling. In addition, the 
larger spin value of Dy {F=J=8 and F=21/2 for bosonic 
and fermionic isotopes, respectively) versus Rb {F=l) 
enhances the magnitude of the Berry's phase — hence 
the strength of the synthetic magnetic field — making the 
quantum Hall regime much more accessible. 

Another important difference between lanthanides and 
alkalis lies within the effective single particle hamiltonian. 
While alkalis' consist only of vector terms in spin space, 
lanthanides' possess sizable tensor terms that generate 
nematic order. By varying the relative strength between 
vector and tensor terms — achievable via Raman and Zee- 
man detuning — one can induce a discontinuous transition 
between spinor states in the (dressed) ground state, an 
effect that can easily be observed. 

For a spin-1 Bose gas, SO-coupling can lead to a stripe 
phase which is a superposition of two spinor condensates 
with different momenta [H |T3l |14] . For large-spin atoms 
such as Dy, the ground state can be a superposition of 
several condensates with commensurate momenta, result- 
ing in a spinor condensate periodic along the direction of 
SO-coupling. The symmetries of these states are conve- 
niently described in Majorana representation as a set of 
2F points on a unit sphere ^20J-[25j, each tracing a differ- 
ent trajectory loop over a period in space — a structure 
we refer to as a Majorana spinor helix in analogy to the 
persistent spin-^ helix observed in condensed matter sys- 
tems [27^. Changes in SO-coupling can also induce tran- 
sitions between Majorana helices of different symmetry. 

The effective hamiltonian of the Raman process: As 
in Refs. [H [S] ITO, we consider two-photon Raman cou- 
pling between different magnetic sub-levels of a spin-/ 
ground state. Coupling is via an excited state with spin- 

induced by an electric field S = 'E-\-'E^ in the presence 
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of a magnetic field B along z, as illustrated in Fig. [TJa), 
where E = ^^xe^^^-^-'^-*) + g^ze'^-^-y-^-^\ The first 
term of E drives cr+ and cr_ transitions with respect to 
the z spin quantization axis, while the second term drives 
the TT transition. 

To derive the effective hamiltonian for the /-spins un- 
der the Raman coupling, we extend the calculation in 
Ref. [28 for purely light to the configuration at hand. 
Eliminating the excited states, we obtain (in the same 
notation as Ref. [28 ) H^i = "^Zf ^ff / ^ff'^^ where 
OLjji^ PfdPf'd^Pf, d is the dipole operator, the P's are 
projectors onto spin / and A/// = u^^a- — ^ff is the 
detuning from the energy difference uoff between the / 
and f levels in a tt or <j transition, and uj^^^cj is either 
uJt^ or ujcr depending on the transition. We decompose H 
into scalar, vector, and tensor components (denoted as 
^R(o), ^R(i), and i^R(2) respectively), 

(0) 2 

1 a\r,g^ - ( \ t t \ 

. (1) 2 

^R(2) = -^^;^{F+(F, + ^)(aJa+A + aLaoA*) 

+F^aLa+ +h.c.}, (1) 

where A = e^[ika+k^)x-\-Auj-Lt] ^ ^ photon creation 

operators for modes t(x ± iy) and z respectively, F± = 
^y^z are spin-/ operators, // is the 
identity matrix, and Aujl = cj^^— cj^r. We have taken g^^ = 
ga = 9 for simplicity. a^fj,{i = 0,1,2) are polarization 
constants given in Ref. [28 corresponding to a rank-z 
{i = 0, 1,2) Raman coupling in the ground state spin-/ 
manifold through an intermediate excited spin-/' state. 

In deriving Eq. ([T]), we have used two conditions: (i) 
the detuning A/j/, is much larger than the difference 
in Zeeman energy both among spin-/ ground states and 
among spin-/' excited states, and thus one can ignore 
the Zeeman-shift-induced m dependence in the detuning 
Afff [29j; and (ii) Aujl <^ uJa^^Tv^ and thus, one can set 
oj^^^ = cJ, the mean value. Both conditions (i) and (ii) 
can be easily satisfied. 

In addition, we assume (iii) that Aul is close to the 
Zeeman energy i.e., Aul ~ cj^, or (5 = Aul — uOz ^ 
ujz' This condition, which is easily satisfied, allows us 
to simplify 11^02) by ignoring the term: the rotat- 
ing wave approximation transforms away the time depen- 
dence of A by performing a rotation g-^^^i^^^^. The 
term gains a phase factor e"^*^'^^^, which then averages 
to zero. 

Comparison between alkalis and open-shell lan- 
thanides: For alkalis, the Raman transitions operate on 




FIG. 1: (a) Raman laser and quantization B-field configura- 
tion; (b) Raman coupling scheme for an alkali atom like Rb 
(left) and for Dy (right); (c) Raman transition energy level 
diagram for magnetic sub-levels of the F = 8 ground state of 
bosonic Dy. 



the ground state and excited states '^P3/2 and 

^Fi/2 possessing a fine-structure splitting Aps- When 
the detuning is much larger than fine-structure splitting 
Aff^ > Afs, ^R(i) is of order g^A^s/A}f, m EQ], 
where A/// = uJ — Fps and Fps is the energy differ- 
ence between the P state and the 5* state without fine- 
structure splitting. The reason is that the states ^P3/2 
and ^Fi/2 couple to the ground state through 
a'^^^ with Clebsch-Gordon coefficients of opposite sign, 
and the /'-sum in ^r(i) in Eq. ([l]) is then of the form 

In contrast, in lanthanides such as Dy, the fine-structure 
splitting of excited state is very large. For instance, 
the 741-nm Dy Raman transition couples to a single ex- 
cited state ^Kg [31 , resulting in a more favorable scaling 

Moreover, Hr(^2) will vanish for the alkalis once Afff 
exceeds the hyperfine splitting. This is because the 
ground state is J = 1/2, which has no matrix ele- 
ment through the rank-2 operator a^^^ back to the same 
manifold. The hyperfine interaction, however, turns the 
ground state of ^^Rb bosons into a spin F = 1 particle, 
leading to H^^2)' ^ ^^Arf/A^^, <C H^^i^'- In contrast, 
the ground state of Dy bosons and fermions have large 
orbital and spin angular momentum, and the matrix ele- 
ment a^^^ is non-zero within the ground state manifold, 
even without a hyperfine interaction: H^^2) ^ d'^/^ff ^ 

^R(l)- 

The heating rates for both alkali and lanthanides 
are given by F ~ g'^Y / ^"ff^^ where 7* is the ex- 
cited state linewidth of atomic species i. Since the 
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FIG. 2: (a) The dispersion of the lowest branch for the Dy 
case with different S [indicated by arrows in (b)]. (b) /cmin 
as a function of 6 for Dy. The dashed hne in (b) excludes 
the tensor term, (c) Majorana spinor representation for two 
ground-state spinor wavefunctions at different 6 [marked by 
(1) and (2)] in (b). (d) dkmin/dS (which is proportional to 
Befi) for Dy (black solid line; black dashed line excludes tensor 
term contribution) and for Rb (orange solid line). In all cases, 
Q is fixed at 4^l, and for Rb case the strength of quadratic 
Zeeman term UqFz is taken as Uq — 1.9^l as used in Ref. 
[1]. In all plots, [/cmin, ^k, 5] are in units of [2A:l, ^l, ^l]. 

linewidth of the ^Kg excited state of Dy at 741 nm 
is ~10~^x narrower than that of alkali's relevant ex- 
cited states, there is much less heating for Dy for the 
same amount of detuning and laser intensity. On the 
other hand, V^^ - {^^^ / /\f f>)H^l^y and r^^^^^^ - 

j^^alkali/^alkali^^alka^^ rpj^^^^ ^-^^ ^^^^ Strength 

of Raman coupling, the ratio of two heating rates 
(pDy/paikaii) _ ( AI^^^V A/// ) (7^^ /7^i^^iO cau be several 
orders of magnitudes smaller than unity for practicable 
laser intensities. 

The single particle hamiltonian for bosonic Dy. With 
the photon fields replaced by their mean values {a±^o), 
Eq. ^ can be simplified to 

Hr = n{e-'^''^yA + h.c.) + 5F,, (2) 
A = F+[if-2C{F, + if/2)], (3) 

where C is a constant that depends on f and /. For 
the 741-nm transition to the ^Kg state of a ^^^Dy atom, 
we have = / + 1. Using the expressions for 
(i = 1,2) given in Ref. [28] and replacing a±^o with their 

expectation value, we find Q = ^ 2A~\^^°^ 
where /3 = (/ + 1) and C = 1/(2/ + 3). 

The term 6Fz arises from the detuning of the frequency 
difference Aujl with the Zeeman frequency ujz- Note, the 
tensor quadratic term F+F^ is absent in previous studies 
of alkali atoms [32l. To understand the physics of Hr^ 



it is useful to apply a spin rotation along z to remove 
the phase factor in Eq. Q [5l [13]. The single particle 
hamiltonian along y is then 

F=(^(^^^^^ + ^F.j+l^(A + At). (4) 

Abelian regime and synthetic magnetic field. The 
abelian regime corresponds to large Q and 5^ and a 
synthetic magnetic emerges when there is field gradient 
5 = Gx pi, n!F . In this regime, Eq. ([2| at each point 
r = (x^y^z) has a unique minimum |'^(r)) in spin space 
separated from other excited states by an energy gap of 
order Q. Projecting H onto this lowest state, one obtains 
the hamiltonian resembling that of a charged particle in a 
magnetic field, where the magnetic field is given by Vx A, 
with the Berry's phase connection A = — i(^(r)|V|^(r)) 
[33] . A quick way to obtain the strength of the synthetic 
magnetic field near r = is to calculate the spectrum 
of Eq. (|4| in semi-classical limit, i.e., by calculating the 
spectrum for a given S and expanding about its minimum 
(0, /cmin((^),0) along ky: 

nky^^) = ^ih-kmUS)f. (5) 

The minimum kmin can be regarded as Ay. The corre- 
sponding synthetic magnetic field is Bq^^ = dAy/dx = 
{dAy/d6){d5/dx) ex dk^in/dS [HIIlj. 

It is important to note that in the abelian limit, the 
ground state of H is controlled by two parameters Vt and 
(5, which drive the system towards the minima of A and 
5Fz^ respectively. As 5 increases, the ground state will 
change from one minimum to another. This can be seen 
in the minimum /cmin(^^) of the spectrum S{ky^S), which 
undergoes a jump at Sc = —6.9 as shown in Fig.[2|a) and 
(b). The nature of the ground state is shown in Fig.[2|c) 
in Majorana representation. For 6 < Sc^ the system is 
ferromagnetic like, with all 2F = 16 Majorana points 
close to the north pole [see Fig. 2(cl)]. (A complete col- 
lapse onto the north pole results in a full ferromagnet.) 
As S exceeds the state gains more nematicity, as seen 
from the migration of Majorana points towards the south 
pole [see Fig. 2(c2)]. (Maximum nematicity corresponds 
to equal distribution of Majorana points at both poles.) 
The dashed line in Fig. [2|b) shows the behavior of kmin 
without the tensor term. In this case, there is no discreet 
jump in /cmin, which is similar to the case of ^''Rb ^16j . 

In Fig. |2|d) we compare dkmin/d6 (i.e., the 5eff 
achieved for the same Zeeman field gradient) between 
Dy and Rb. With the tensor term excluded (dashed 
line), the maximum value of dkmin/dS achieved for Dy 
is larger than that for Rb case by a factor of ^F = 8. 
This is because the larger spin structure enhances the 
maximum momentum transfer from AkL in Rb to 4F/cl 
in Dy. Including the tensor term (black solid line), the 
synthetic magnetic field dkmin/d5 is greatly suppressed 
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when 6 < 6c, but is much enhanced for 6 just above 6c. 
Operating in the latter range of 6 wih provide a very 
strong synthetic magnetic field. 

SO- coupled regime and Majorana spinor helix: The 
single-particle dispersion exhibits multiple minima that 
are almost degenerate in the SO-coupled regime. This oc- 
curs for small 6 and Q. For ^1 = 0, Eq. Q exhibits 2F + 1 
degenerate minima with wavefunctions e^^'^^^^ Lp^ where 



For smah the 2F 



im2ki^y ^ 

1 local minima 



remain, but their energies are no longer degenerate 
and their wavefunctions are modified to (in leading order 

of n): 



l^m^m+l) 



(6) 



where /^^ = Vt{m ^ l\K\m) . In the presence of interac- 
tions, the ground state will be a linear combination these 
modified minima <l> = ^^(.mi^m- Expressing <l> in the 
original basis (/p^, which are spin components measured 
in a Stern- Gerlach experiment, we have ^ = X^^^m^m^ 
where 



Jm2kLy 



i-e 



-i2kLy , 



(7) 

The coefficients ^rn can be determined in a straight- 
forward fashion by minimizing the energy 8 = 

Z)m=-F^^I^^P + ^int({^m}), whcrC fint({^m}) is the 

interaction energy of the form 

j=0 m=-j 

Y{^J,m\J,m\ J,m - m')r]rn'r]m-m^ , (8) 
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FIG. 3: (a) Single particle spectrum for ^ = and Q = 
0.083^L- Ek are in units of and El, respectively, (b) 
A one-parameter phase diagram in term of Q/Ej,. 6 is fixed 
at zero. For the purpose of illustration, we choose a set of pa- 
rameters {^0, • • • , ^'le} = {6, 7, 7, 7, 7, 2, 0, 7, 7} as an example 
(other interaction parameters will result in spinors of differing 
symmetry, see Ref. [25 ; these parameters remain unmeasured 
for Dy). With this set of interaction parameters, the system 
is in an octahedron Oh phase at zero Q, and Csvia), Czv{fi) 
and Coov phase for < Q./Ei^ < 0.101, 0.101 < Q/Ei^ < 0.133 
and Q/Ei, > 0.133, respectively. The phases are labelled by 
symmetry in the dressed state basis, (c) Schematic of a "Ma- 
jorana spinor helix" in C3v{l3) phase of (b). The grey lines are 
the trajectory of each Majorana point on the sphere versus y 
translation. 



and g2j are j + 1 independent interaction parameters in 
different total spin channels [34] , and rjm is related to ^rn 
via Eq. The ground state is sensitive to variation of Q 
since Raman coupling effectively changes the interaction 
between dressed states; in Fig.[3]^b) we give an example of 
the change of symmetry of the ground state as a function 
of ft [35 . 

Equation ([t]) shows that each spin component of the 
ground state <I> = T]m^m is a periodic function in y 
with wavevector /c^. Such variations can not be described 
by spin rotations. Rather, it corresponds to each of the 
2F Majorana points tracing out different loops as the 
atoms travel over a period ir/kL along y, as shown in 
Fig. 3(c). We name such a structure a "Majorana spinor 
helix," as opposed to a reshuffling of different Majorana 
points after a period, which would be (in the analog of 
superffuid ^He) a Majorana "soliton." 

Concluding remarks: The outstanding challenge for ex- 
perimental research employing a synthetic gauge field is 
the reduction of heating. Our discussion show that this 
problem can be solved by using atoms with large orbital 
and spin angular momentum such as Dy, which also has 



the advantage of considerably increasing the strength of 
the synthetic gauge field. Moreover, SO-coupled high- 
spin bosons and fermions possess far richer classes of bro- 
ken symmetry ground states than spin-1/2 fermions and 
spin-0 and spin-1 bosons. The realization of a strong 
synthetic gauge field and SO-coupling in these systems 
of large-spin atoms will constitute a major step toward 
exploring exotic correlated states of fundamental impor- 
tance that are difficult, if not impossible, to realize in 
solids. 
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